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Abstract 

Particle trapping in multi-state quantum walk on a circle is studied. The time- 
averaged probability distribution of a particle which moves four different lattice sites 
according to four internal states is calculated exactly. In contrast with "Hadamard 
walk" with only two internal states, the particle remains at the initial position with 
high probability. The time-averaged probability of finding the particle decreases 
exponentially as distance from a center of a spike. This implies that the particle 
is trapped in a narrow region. This striking difference is minutely explained from 
difference between degeneracy of eigenvalues of the time-evolution matrices. The 
dependence of the particle distribution on initial conditions is also considered. 
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1 Introduction 



The quantum walks attract our attention not only in a fundamental field of 
quantum mechanics [1] but also in a filed of quantum information [2] (see also 
recent review with many references [3]). This is because the quantum walks 
will be achieved on quantum computers [4], and it may be used as a new 
quantum algorithm. 
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A classical random walker on a line moves to a nearest left site or a nearest 
right site randomly at every time step. On the other hand, a quantum walker 
moves, in superposition, both left and right according to additional degree of 
freedom called "chirahty" . First of all, we describe the Hardamrd walk, which 
is the most well-known quantum walk. The Hadamard walk is characterized by 
two parameters. The first is the spatial position and the second is the chirality, 
which takes values and "—1". The corresponding Hilbcrt spaces are 

H = and Hj = C^, where C is the set of complex numbers. The position 
space consists of basis states \x) G H{—oo < x < oo) corresponding to the 
walker located at n-th lattice site. Possible experimental implementations of 
the Hadamard walk have been proposed by a number of authors [4-6]. For 
example, a quantum walker is an atom and internal states of the atom such 
as hyperfine structure are used to distinguish the chirality. 

Let ip{s, X, t) be a value of wave function at position x and time t with 
the chirality s G { — 1,1}. The wave function is changed by superimposing 
ip{—s, X — s,t — 1) on ip{s, X — s,t — 1) at ever time step (details are defined in 
the next section). Assume that a walker exists at the origin in spatial-temporal 
plane, then we find the walker in the symmetrical Hadamard walk with the 
same probability at the position 1 or —1 after a single transformation. This 
behavior is the same as symmetric classical random walk. The primary ques- 
tion in the quantum walk is to determine the probability distribution of the 
walker after repeating the transformation of wave function. The recent exten- 
sive studies made clear mathematical properties of the Hadamard walk. We 
know that the probability distribution is quite different from that of classical 
random walk. When the initial wave function is given by 1,0,0) = 1/2 
and 0,0) = i/2, (symmetric case), we see that double peaks spread in the 
mutually difi'crcnt direction. Furthermore the standard deviation increases in 
proportion to time. This means that the quantum walker spreads more quickly 
than classical random walker, and it is a desirable feature for seeking problems. 

Recent studies revealed that there are two distinguishing profiles in modified 
Hadamard walks on an infinite space. When the walker exists initially at the 
origin, the one of shapes is concave in a broad way and the other shape is a 
spike. The spike at the initial time in the symmetrical Hadamard walk sprits 
two peaks and they move away from the origin as time increases. The proba- 
bility of finding the walker near the origin decreases and approaches to zero. 
Thus the probability distribution of the original Hadamard walk belongs to the 
former. On the other hand, we can always find a walker with high probability 
near the origin in the later case. We know little about the spiky probability 
distribution in the later case. The existence of the spikes are already reported 
in the two-dimensional Grover walk [7,8], one-dimensional quantum walk as 
many coins [9], but, no probability distributions are calculated explicitly. In 
this paper we consider quantum walk with four internal states, which is the 
same as the quantum walk with many coins essentially, and compute analyti- 
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cally the probability distribution. The Hadamard walker described above has 
two internal quantum states and moves to only nearest neighbor sites, while a 
walker in "four-state quantum walk" can move beyond nearest neighbor sites 
by a single transformation. 

The rest of this paper is organized as follows. We define multi-state quantum 
walk generally and show a result obtained by simulation in the next section. 
In section 3, eigenvalues of the time evolution matrix are exactly calculated 
and the wavef unction is expressed using those eigenvalues. In section 4, we 
introduce a time-averaged probability and show that the quantum walker is 
found with high probability at the initial position. In section 5, a quantum 
walker with three internal states, in which a walker jumps asymmetrically, is 
considered and it is shown that localization is not observed. 



2 Multi-state quantum walk 

We consider a quantum walker with n different internal states and moves to 
n different lattice sites according to the internal states. Assume that n is an 
even integer larger than 2. Let s E T = {±1, ±2, • • • , ±n/2}. We then define 
the transformation of n-state quantum walk by 



where Sg^s'^s = 2/n and 6s^s = 2/n — 1. This transformation is quite simple, 
but, it produces rich quantum effects. If the number of internal state is four, 
we regard the above transformation as the quantum version of the classical 
random walk in which a walker jumps to each nearest site and next nearest 
site with the same probability 1/4. 

Our main conclusion arguing that the particle in n-state quantum walk is 
trapped near the starting point is available for any even n except n — 2, but 
in this paper we focus on four-state quantum walk to avoid excessive generality. 
The transformation with n = 2 is diffrent from the Hadamard transformation, 
and the dynamics of 2-state quantum walk is trivial. 

Wc begin to calculate the probability distribution of the quantum walks on a 
circle containing N sites [10]. We define the probability of finding a walker at 
position X and time t by 



l/j{s, X,t + 1) = ^s,s'1p{s', X 



s,t), 




s'er 



P{x,t) = J2ij{s,x,t){i;{s,x,t)y, 



(2) 
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where "*" denotes complex conjugate. To compare the probabihty distribution 
between the Hadamard walk and four-state quantum walk, Fig.l shows the 
snapshot of probability distribution at t = 50 starting from an initial wave 
function: 



ip{s, x,0) — < 



1 

a/" 
i 





X — ^ and s < 0, 
X — ^ and s > 0, 



(3) 
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Fig. 1. Comparison of probability distributions between the Hadamard walk and 
four-state quantum walk at t = 50. A walker initially exists at the origin in both 
cases. The walker in the Hadamard walk moves away form the origin. The walker 
in four-state quantum walk, by contrast, stays at the origin with high probability. 

In contrast with the Hadamard walk we clearly find a single spike at the origin 
in four-state quantum walk. A similar spike has also been observed in two- 
dimensional quantum walk by simulation [7] and the height of time-averaged 
probability at the origin is exactly calculated [8]. However, the probability 
distribution has not been obtained explicitly as a function of the location. For 
these reason, the main purpose of this paper is to calculate the time-averaged 
probability distribution of four-state quantum walk and show that the spike 
in Fig.l does not disappear even in the thermodynamic limit. 



3 Transfer matrix and eigenvalues 



The time evolution of the wavefunction of the multi-state quantum walk ^(t) is 
determined by transformation Eq. (1) and it is expressed briefly by introducing 
a time evolution operator M satisfying 1) = M^'(t). Since the matrix M 
is unitary, '^{f) is explicitly given by transforming M to a diagonal matrix as a 
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function of the time. Using Fourier transformation it is found that eigenvalues 
of M are calculated from the eigenvalues of the following matrix: 



u-'^^ 

0;-^= 

cu'^ 

a;^'^ 



-1111 

I- 111 

II- 11 

III- 1 



(4) 



where uj = e^'^^^^ . Changing the value of k from to iV — 1, we get all 
eigenvalues of M (see details in Ref. [11]). The eigenvalues of are given by 



Afe,i — —1, 



Afc,2 — 1, 



- _-9k-i^J'^- gl - _ -9k + i^J'^- gl 

^k,3 — 7^ , '^M — o • 



(5) 
(6) 



where = cos 62k + cos^4fc and O^. = krr/N for A; = 0, 1,---,A^ — 1. We stress 
here that the eigenvalues —1 and 1 are independent of value k, put simply, 
they are strongly degenerative. In fact, the degree of degeneracy for —1 and 
1 are N + 2 and N, respectively. Since Qk — QN-k for k > 0, the degree of 
degeneracy of other eigenvalues is two unless N is a multiple of 5. In case 
of = even, the classification of eigenvalues is more complicated. Thus we 
suppose that N is odd and is not a multiple of 5 in this paper. 

All eigenvalues are now obtained and then the wave function is formally ex- 
pressed by 



X, t) = ci;i,(x)(-i)* + c'j'^x) 

(N-l)/2 {N-l)/2 

+ E c,,,,3,jv(:^:)e^"'='^* + E c,,fe,4,iv(^)e^"'=■*^ (7) 



J2) 



fe=l 



fc=l 



where u!k,n denotes the argument of Xk,n, and the first and second terms express 
the contribution of eigenvalue —1 and 1 to the wave function. The coefficients 
in the above formula are calculated from the eigenvectors of M. Let (f)k,n be 
the eigenvectors corresponding to the eigenvalue „ of and (f)k,n.j be the 
j'-th element for i = 0,1,2,3. We assume additionally ■ ■ ■ , 0^,4} is an 

orthonormal basis. Then the i-th element of eigenvectors of the matrix M 
with Afc,n, which becomes an orthonormal basis is obtained by 



^k,n,j - 



■'k,n, J mod 4'-^ 1 



(8) 
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where [^J denotes the integer part of z. If the walker exists initially at the 
origin, the coefficient Cs,hfi,N{^^i for example, is expressed by 

1 ^ 

Cs,fc,3,Jv(2:) = — ^{(^'^"0fc,3,.(0Mj)* + ^^'^"''^"'/'iV-M,3(0iV-fe,3j)*}V'(s,-,O,O), 

(9) 

where s\ = —2, S2 = —1, S3 = 1, and S4 = 2. 



4 Time-averaged probability 

In general, the value of wave function at any fixed position x does not converge 
to a constant value. Thus we here define the time-averaged probability of 
finding the particle at the position x and in the internal state s in order to 
show the localization in long time scale as 



1 

Pjv(s,a;)= lim -J2^Nis,x,t){'iljN{s,x,t))*. (10) 

Noting that hmr^oo E^V e^^'^-'^'^VT = for a; ^ uj' , we have the following 
equation: 



Pn{s,x)^ c%{x)^ + c^^^l{x)^ + {\cs,k,z,N{^)? + \(^s,kAN{x)\^)- (H) 



,(2) 



{N-l)/2 



k=l 



Let us consider time-averaged probability in the limit N ^ 00 below. Com- 
bining Eqs. (8) and (9), we find that |cs,fc,3,Ar| and |cs,fe,4,Af| decrease in the form 
1/N for large N. Therefore the third and fourth terms vanish in the limit of 
iV — cxD. As a result, time-averaged probability on a circle containing infinite 
sites is simply given by 



Poo{s,x) = lim 

JV— >oo 



x) 



+ 



J2) 

''s,N 



[x) 



(12) 



Assume that the particle exists at the origin initially. Then the values of 
c^P^{x) for / = 1,2 in the right-hand side are determined from only four 
components of the initial state and they are expressed in the following form: 
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—z,oo\ / 




C-l,oo(^) 




^1,00 (■^) 









U%.,{x) U%_,{x) U%,{x) U%,{x) 

U%_,ix) U^L.ix) U^ix) U%ix) 



^(-2,0,0) 
^(-1,0,0) 

VXi,o,o) 

^^(2,0,0) 



(13) 



In this formula, Uf\ expresses a contribution from the initial state '?/'(sj, 0, 0) 
to c'^P^ix). If each component of the above matrix U^^"^ in Eq.(13) is exactly 
obtained, then the time-averaged probability is immediately calculated from 
Eq.(12) for any initial state. 



Let us now calculate UgJ{x). Using the set of eigenvectors of iJfe, Ugj^{x) is 
given by 



(0, 



Ufl{x)^ lim ^ V 



N-l 



(14) 



fe=0 



We use the following orthonormal eigenvectors corresponding to \k,i for / = 1, 2 
and A; > 0: 



0, 



k,l 



m) 



(15) 



where 



/i(^) = V6T4cos2^, 

f2{e) = 1/2(1 + cos 2^) (2 + cos 2^ + cos 4^), 

= -6,4(^fc) = 1, 
ei,2(^fc) = -6,3(^fc) =2c0S^fc, 
6,l(^fe)=6,4(^fc) =1+C0s(2^fc), 

6,2(^fe)=6,3(^fc) =cosefe + cos(3efe). (16) 



In case of /c = 0, since eigenvalues are {—1,1,— 1,-1}, and the value — 1 
degenerate, we use the following eigenvectors which are orthonormal basis 
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00,2^ 
00,3^ 



1 1 \ 

1111 

2' 2' 2' 2 

^/6'^/6' ' VqJ' 

1 1 Vs 

~27^'~275' 2~' 



2x/3> 



(17) 



Plugging above eigenvectors into Eq. (8) and carrying out the summation 
in Eq. (14) we have UgJ(x) exactly after a somewhat tedious calculation as 
follows: 



Ui^hx)^g?i\x\), 



(0/ \ 
92 [x) 



where 



9?{l-x) 



9f 



x > 1, 

X <l, 

X = 0, 
otherwise. 



(18) 



9f\x) 
9?\x) 



a 



cos(7a; + 5i) 
2l7^ ' 



20 



5 + ^ (1) 5-\/5 

-a , q\ ' = — , 

' 10 ' 

cos(7a; + ^2) 



9?\x)^ 



1 V14 + 28\/2 

2 28 ■ 



The value of a, /?, 7, 5i and 82 are given by 



(19) 



« = ^""2"7! + ^ 2"^' 7 = tan-^(Vll + 8v^ 



tan ^ "\/7 
oi = -z , 02 



( 



tan 



V 



3 + 2 



^ 2^206 + 148^2 



(20) 



When the number of size is large, the valuable 6'^ = ki^/N is almost con- 
tinuous. Thus the infinite summation is calculated by transforming the sum- 
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mation into the complex integration on the unit circle. The integration is 
exactly calculated by residue theorem, and a and (5 are poles of the function 
of u!''^4>k,i,s{4>k,i,i)* in Eq. (14), in which 9k is replaced as continuous valuable. 
Other components are obtained from the following relations: 

f/i5(,T) = + i^-1,-2, 

uS{x)^{-iyU%{x-s) 5 = -1,-2. (21) 

Although the obtained expression is rather complicated, this formula gives 
exact values of the time-averaged probability for any initial state and for any 
position of lattice site. 

Let us write briefly the initial state by — ('0(~2, 0, 0), • • • , ■?/'(2, 0, 0)), 
and consider the time-averaged probability summing over all possible states 
Poo{x) = J2ser Poo{s,x). Fig.2 shows Poo{x) for three different initial states: 
(a) symmetrical case ^'o = (|, |, |, |), (b) asymmetrical case \l'o = (1, 0, 0, 0), 
and (c) extinction case = {—^, —^). One clearly finds that the walker 
is localized near the origin for the case (a) and its shape is almost the same 
with that in Fig.l. In case of (b), there is a double peak and it shifts to the left. 
In contrast with (a) and (b), no peaks exist in the case of (c). The probability 
of finding a walker decreases in inverse proportion to the system size and it 
converges to zero in the limit of A'^ — > oo. 




-4-2 2 



X 

Fig. 2. Dependence of the time-averaged probability distribution of the four-state 
quantum walk on the initial states "^q: (a) (^,^,|,|), (b) (1,0,0,0), and (c) (-^, 

2'2' 2J- 

We move our attention to the asymptotic behavior for large \x\. From Eq. (19) 
each component of the matrices C/^'^ exponentially decreases as increasing the 
\x\. The values of a and f3 are nearly equal to —0.381966 and —0.53101, re- 
spectively, thus the latter contributes largely to Poo{x). Fig.3 shows logg -Poo (2;) 
as function of x with the initial state (a). The slope of inserted line means 
the exponent of decay is —1.266. This exponential decay of the time-averaged 
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Fig. 3. Semi-log plot of the time-averaged probability of four-state quantum 
walk with initial state (a). The dashed straight line in the figure shows that the 
time-averaged probability decays exponentially for large x. 

probability Poo{x) shows surely that four-state quantum walker localizes near 
the origin except for special initial states. Very recently similar exponential de- 
cay is observed in momentum space of generalized quantum walk by Romanelli 
et al. [12]. We should remark that only the four-state quantum walk starting 
from ^'o = {— e*^, e*"^, e*^, — e*'^}/2 for any real number Lp does not show the lo- 
calization. This means that the localization is substantially universal property 
in four-state quantum walk. 



5 Three-state quantum walk 

We have shown that the four-state quantum walker is fixed near the initial 
position except special initial conditions. In this section, three-state quantum 
walk with asymmetrical jump is considered. The time evolution of three-state 
quantum walk is defined by 

^(-l,a;,t + l) = ^{-^(-l,a; + l,t) + 2^/'(l,a; + l,t) + 2^/'(2,a; + l,t)}, 

^'(l, X, i + 1) = ^ {2V(-1, x-l,t)- ^^(1, X - 1, + 2V'(2, x-1, t)} , 

^(2, x,t + l) = ^ {2^(-l, x-2,t) + 2^(1, x-2,t)- ^(2, x-2,t)}. 

(22) 



In four-state quantum walk, the walker can jump symmetrically to the right 
and the left. This three-state quantum walker, by contrast, can jump to nearest 
neighbor sites and a right-hand next-nearest neighbor site, but the walker can 
not jump to a left-hand next-nearest neighbor site. 
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Fig.4 shows snapshots of probabihty distribution at t = 100, 200 and 300 
starting from the origin with an initial state ■?/'(— 2,0,0) = |, tjj{^l,0,0) = 
0, 0) = |. We observe that the several local peaks move to the right and 
the heights gradually decrease in the prolonged simulation. 
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Fig. 4. Probability distributions of three-state quantum walk at t = 100, 200 and 
300. A spike, which exits initially at the origin, is divided into several packets and 
each packet moves to the right. The height of each peak becomes low as time passes 
and approaches to zero in the limit oi t ^ oo. 

The striking difference between four-state and three-state quantum walks is 
explained from difference in the matrix if/j. The matrix corresponding to 
three-state quantum walk is given by 
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uj-^ 







-1 


2 


2 


uj^ 







2 


-1 


2 









2 


2 


-1 



(23) 



The eigenvalues of h)^^ are exactly obtained. It is easily found that the each 

degree of degeneracy is independent of the system size. In other words, there 
is not any eigenvalue corresponding —1 and 1 considered in the previous sec- 
tion. Therefore the height of each peak inevitably approach as to zero in the 
thermodynamic limit. Consequently, the localization can not be observed in 
the three-state quantum walk. 



6 Summciry 



The key factor of the localization in the four-state quantum walk is the de- 
generacy of eigenvalues. If the degree of degeneracy is independent of the 
dimension of the Hilbert space, the contribution to the time-averaged prob- 
ability approaches to zero in the limit of oo. Therefore the existence 
of eigenvalues whose degree of degeneracy is proportion to the size of the 
Hilbert space is a necessary condition. This basic condition is also satisfied 
in the two-dimensional Grover walk. As a result, the walk shows localization 
[8]. Prom this criterion we conclude that the Hadamard walk does not ex- 
hibit localization. Furthermore we proved that three-state quantum walk with 
asymmetrical jump does not show the localization. 

The time-averaged probability at the initial position in the Grover walk al- 
ready calculated. However the probability distribution has not been calculated 
yet. Thus we showed firstly that the time-averaged probability of finding the 

particle decreases exponentially (not Gaussian) as distance increases from a 
center of a spike and there are special initial states with which the spike dis- 
appears in a one-dimensional quantum walk. 

The primary criterion of being localization in probability distribution is only 
the dependence of degeneracy of eigenvalue on the system size. If the number 
of the internal state is even, then eigenvalues include —1 and 1 and their 
degree of degeneracy increase as the system size. Accordingly we can conclude 
that the multi-state quantum walker with even internal states almost freezes 
without detailed calculations. 
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